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We investigate the relation between the existence of mass inflation and model parameters of string- 
inspired gravity models. In order to cover various models, we investigate a Brans-Dicke theory that 
is coupled to a (7(1) gauge field. By tuning a model parameter that decides the coupling between 
the Brans-Dicke field and the electromagnetic field, we can make both of models such that the 
Brans-Dicke field is biased toward strong or weak coupling directions after gravitational collapses. 
We observe that as long as the Brans-Dicke field is biased toward any (strong or weak) directions, 
there is no Cauchy horizon and no mass inflation. Therefore, we conclude that to induce a Cauchy 
horizon and mass inflation inside a charged black hole, either there is no bias of the Brans-Dicke 
field as well as no Brans-Dicke hair outside the horizon or such a biased Brans-Dicke field should be 
well trapped and controlled by a potential. 
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I. INTRODUCTION 


Motivated from quantum gravity m as well as cosmology [1] and black hole physics [3], many models 
of modified gravity have been suggested. One big branch of the modified gravity is the scalar-tensor gravity 
013 that can be motivated from string theory [8 . In these classes of string-inspired models, the investigation 
of black hole physics is very worthwhile, especially not only for static solutions [9|, but also for dynamical 
behaviors. 

For this purpose, in the authors previous paper m, we investigated the prototype of the string-inspired 
model as 


s= wih'^ d,x 


$ 1 ? - 


( 1 ) 


where R is the Ricci scalar, $ is the Brans-Dicke field that can be interpreted as a dilaton field (if u> = — 1, 
while there can be other string-inspired values for w El), and F is the two-form field that can couple to a 
complex scalar field. In order to investigate dynamical properties, we implemented a numerical formalism, the 
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double-null formalism m, and people have obtained various results mm. Especially, in the authors’ first 
paper, we mainly focused to see the causal structures and responses of the Brans-Dicke field <f>. 

One important result of the authors’ previous paper was to related causal structures of string-inspired charged 
black holes (see also m) with the responses of the Brans-Dicke field mm Because of a certain coupling 
between the gauge field and the Brans-Dicke field, the responses of the Brans-Dicke field during and after a 
gravitational collapse can be sensitively changed. Because of this, the internal causal structure of a charged 
black hole could be sensitively depend on the choice of coupling parameters ui and /3, whether there exists a 
Cauchy horizon and mass inflation |18l or not. In conclusion, we observed that if w > —3/2 and /? > 0, then 
there was no Cauchy horizon, where this result is consistent with m- 

However, the next natural question is this: is that all that determines the internal structures? In fact, 
there are much more complications in realistic string-inspired models. For example, let us choose ui = —1 
to discuss a dilaton model. However, to be consistent with the experimental results [20], such a dilaton field 
should be stabilized by a potential, so to speak H(d>). Then this potential will restrict the responses of the 
Brans-Dicke field. Then will it change the internal structure of a charged black hole so that there appears a 
Cauchy horizon again? Like this, we can ask whether a certain choice of a string-inspired model (e.g., choosing 
/3) solely determines the internal structure, or there are lots of choices of parameters that determine the internal 
structures. This paper will be devoted to answer on this problem. 

This paper is organized as follows. In SEC. |TIJ we discuss a detailed model of the Brans-Dicke theory with a 


U{ 1) gauge field, where we further introduce a potential term of the Brans-Dicke held. In SEC. Ill we investigate 
numerical results for the existence of mass inflation and Cauchy horizons by varying some parameters, especially 
the coupling parameter /? and the mass scale of the potential. In SEC. |IV[ we summarize our results and discuss 
on possible future issues. In Appendices A, B, and C, we discussed details of our numerical formalism, initial 
settings, convergence and consistency checks, and some details of numerical results. 


II. MODEL FOR CHARGED BLACK HOLES 

In this section, we discuss on the details of the model for a Brans-Dicke theory with a t/(l) gauge held. In 
addition, we briefly summarize the previous results of the authors m and discuss on the motivation of this 
paper. 


A. Brans-Dicke theory with a U( 1) gauge field 

The prototype action of the Brans-Dicke theory with a U( 1) gauge held becomes (c = G = h = 1) 

1 


^ = / V^gd 4 


16tt ( $i? “ - V($)) + ^£ EM 


( 2 ) 


where 


£ EM = -- (</ ;M + ieA^)g^ (0 ; „ - ieA v 4>) - —F^F^ 


( 3 ) 
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FIG. 1: Conceptual interpretation of the results of double-null simulations. Left is the numerical integration domain. 
By tilting 45 degrees, we can interpret this as a Penrose diagram. 


4> is the Brans-Dicke field, u> and /3 are free parameters that determine the model, <f> is a complex scalar field 
with a gauge coupling e, A M is a gauge field, = A v .p — A M; „, and is the potential of the Brans-Dicke 
field, where we use the form 

D(4>) = ^M 2 ($ - l) 2 (4) 

with a constant M. 

The Einstein equation becomes as follows: 

= 87t (T™ + &- l TV v ) , (5) 


where the Brans-Dicke part and the matter part of the energy-momentum tensors are 


rriBD _ 1 

^ ~ 8n<f> 


( 9nv&;pag P + ^;p,iy) + g 7r( j,2 ( 0 9 pi* , p*&\o g P ) 


167r$ 


9pw> 


T pu = 9 + X (-<!>■,PieA v (j) + < p.JeA^ + <j>.pieA„</> - ^.JeA^) 


+ - 7 ^FppF v p + e 2 A^A^cfxj) + C EM g^ u . 
47r 


In addition, the field equations for matter fields are as follows: 


8tt4>' 3 


1 


0 = (T c - 2/?£ EM ) - ($D'($) - 2V ($)), 


3 + 2 uj 


3 + 2c j 


p. 


o = FicA^ (2cj). l _ l + ieA^(l)) + ieA^ v g^ v (j)+ — (0 ;i , + ieA„(f>) g^ v , 

0 = ^ (f v p- v + “ ie( t> (</v - (<j)-p + ieAp(j >), 


where T c = T c, ) ; . 

h 1 


( 6 ) 

(7) 


( 8 ) 

(9) 

( 10 ) 
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In the double-null formalism, we use the double-null coordinates 

ds 2 = —a 2 (u, v)dudv A r 2 (u, v)d£l 2 (11) 

and present every equations using this metric ansatz. Here, we assume the spherical symmetry, u is the retarded 
time, v is the advanced time, dVl 2 = d6 2 + sin 2 9dip 2 , where 9 and p are angular coordinates. The detailed 
formulations and assignments of initial conditions are discussed in Appendices A and B. We need to assign the 
initial conditions for the in-going and out-going null slices. Especially we prepare a condition for a collapsing 
shell along the initial out-going null slice. After we get a two dimensional data by solving numerical simulations 
as functions of u and v, by tilting 45 degrees, we can interpret them as a Penrose diagram (FIG. [l]). 


B. Summary of previous results 


In the previous paper of the authors m, we investigated causal structures and responses of the Brans-Dicke 
field of string-inspired models. One of the crucial point is the dynamics of the Brans-Dicke field that satisfies 
(M = 0) 


V 2 $ 


3 -f 2 oj 3 A 2uj 


x Q = 0, 


where 


K. = 8n<f> p 


wz A zw A iea{zs 

2'KOt 2 


Q = 8t 


47t r 4 



( 12 ) 

(13) 

(14) 


Here, Q can be interpreted as a charge term that remains after gravitational collapse, while K. is a kinetic 
term of the matter field that only contribute during a gravitational collapse and disappear as time goes on. 
Therefore, if we assume oj > —3/2, then 


- During gravitational collapse, via the K. term, 

if /3 > 1, then the Brans-Dicke field towards a weak coupling limit, 
if p < 1, then the Brans-Dicke field towards a strong coupling limit, 

- After gravitational collapse, via the Q term, 


if /3 > 0, then the Brans-Dicke field towards a weak coupling limit, 
if /3 < 0, then the Brans-Dicke field towards a strong coupling limit. 


Therefore, for string-inspired models with 0 < /? < 1, after gravitational collapse, the Brans-Dicke field towards 
the weak coupling limit; and, in these cases, as we will repeat again in the next section, there is no Cauchy 
horizon. On the other hand, if there is a tendency that the Brans-Dicke field is not biased (/3 = 0), there exists 
a Cauchy horizon. 
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C. Task of this paper 


Then can we conclude that, for Type I (/? = 1/2) or Heterotic model (/? = 1), always there is no Cauchy 
horizon? However, in general the situations are complicated. In this paper, we vary more parameters. 

— /3 dependence: Even though there is no good corresponding model in string theory, as a theoretical con¬ 
sideration, we can investigate for /3 < 0 or (3 > 1 limits. For example, if /3 < 0, then during and after 
gravitational collapses, the Brans-Dicke field will be biased to a strong coupling limit. Then can there 
be a Cauchy horizon again? In other words, does the tendency of the bias of the Brans-Dicke field solely 
determine the existence of the Cauchy horizon? 


— Potential dependence: In general the dilaton field should be stabilized by a potential. This gives a cor¬ 
rection term to the equation so that 


V 2 $ + x JC 

3 + 2 lo 


3 + 2 w 


x Q- 


1 


3 + 2lj 


x? = 0, 


(15) 


where 


V = ®V' - 2V (16) 

is a contribution from the potential. For any potential, around the local minimum, it can be well approx¬ 
imated by a quadratic form V = M 2 ($ — l) 2 /2, and hence V = M 2 (*I> — 1). Then what is the relation 
with the existence of a Cauchy horizon and the mass parameter M? 


III. MASS INFLATION AND EXISTENCE OF CAUCHY HORIZONS 


In Einstein gravity, mass inflation is a general phenomena that appears inside a charged black hole (18] . For 
four dimensional cases, a static charged black hole metric is 


ds 2 = - 1 - 


2m 


+ \ ] dt z + ( 1 - 


—+ C) ' dr' 2 + r 2 ci0 2 . 


(17) 


In this case, there are two horizons, one is an outer event horizon and the other is an inner Cauchy horizon. If 
there is a pulse of any matter or energy along the in-going null direction, then it can be well approximated by 
T/;, ~ d^vd u vF(v)/r 2 , where F(v) is the luminosity function that is proportional to v~ p with a constant p for 
astrophysical matters. Now if an observer falls into a black hole and approaches the Cauchy horizon following 
the increasing v direction, then the observer measures the energy density p = T IJV l IJ i v ~ e 2K ~ v , where l^ is a 
tangent vector of the observer and K- is the surface gravity at the inner horizon. Therefore, the local energy 
density should exponentially increases as an observer approaches to the inner Cauchy horizons; hence, as a 
back-reaction, relating curvature quantities should increase exponentially around the inner horizon. 

In general, since the curvature quantities are increased seriously, the back-reaction to the metric should be 
very large and hence we need to rely on numerical calculations to see the detailed internal structures |15j . If 
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CO— 1, p=0, M=0 




FIG. 2: r (left) and $ (right) for /? = 0 (top), 0.5 (middle), 1 (bottom) with M = 0. Thick black curves denote a 
space-like singularity (r = 0), thin black curves denote r = const, contours, red curves denote r,„ = 0 apparent horizons, 
and blue curves denote r tU = 0 horizons. 


there are some contributions from higher curvature corrections (51) or a black geometry has a different topology 
[55] , then the tendency of mass inflation can be changed. 

In this section, we see more details on the relations of the Brans-Dicke Held and the existence of a Cauchy 
horizon. The first question is this: does the biased direction of the Brans-Dicke field determine the existence 
of the Cauchy horizon? 


A. String inspired models with the M = 0 limit 

First, we report on the case for /3 = 0, 0.5,1 with M = 0 (FIG. [2]) where the other parameters are fixed and 
commented in Appendix B. As we commented, if /3 > 0, after the formation of the black hole, the Brans-Dicke 
field is forced toward the weak coupling region (i.e., $ > 1), and in the end, there becomes no inner Cauchy 
horizon. 



























B. Variation of parameters 


1. /3 < 0 limit 

Now let us vary the parameter /3 less than zero. Then as we discussed, we will surely see the bias of the 
Brans-Dicke field toward the strong coupling limit after a formation of a black hole. Then can this guarantee 
the existence of a Cauchy horizon? 

To check this, we study the cases of j3 = —1, —0.5 with M = 0 (FIG. [3]). Although we are interested in the 
case u) = — 1 (dilaton limit), to amplify the dynamics of the Brans-Dicke field, we also investigate the cases of 
u> = —1.4. 

These results are impressive. Let us summarize important points: 

- As we can see the Brans-Dicke field, during and after the gravitational collapse, the Brans-Dicke field is 
biased toward the strong coupling regime. 

- For some cases (first, second, and fourth in FIG. |3j) , inside the event horizon, there exists not only r tV = 0 
horizon but also r u = 0 horizon. This means that the in-going observer will see an increasing areal radius 
(left of FIG. [4]) . This is a kind of wormhole inside the black hole. This is extremely surprising, but we 
can understand this. In the Einstein frame, the spherical area is proportional to oc <I>r 2 , where $ and r 
are calculated by the Jordan frame. Therefore, in this case, as an observer falls into the black hole, the 
observer experiences r —> 00 and $r 2 —> 0 (right of FIG. [4|. Therefore, in the Jordan frame, the thick 
black curves in FIG. [3] may not be r = 0, but we should regard the region as the singularity. 

- If we understand this, then in fact, all of these results shows that there is no Cauchy horizon and mass 
inflation. The causal structure is the same as that of the neutral black hole and the cases of (3 > 0. 

Therefore, we can conclude that the existence of the Cauchy horizon and mass inflation is not determined by 
the direction of the bias of the Brans-Dicke field. After a charged black hole is formed, at once it is biased 
toward strong or weak couping regime, then there is no Cauchy horizon. Therefore, this strongly indicates 
that the existence of a Cauchy horizon and mass inflation is related the fact that there should be no bias of 
the Brans-Dicke field; e.g., there should be no hair of the Brans-Dicke field outside the horizon. Perhaps, the 
existence of a Cauchy horizon and mass inflation is the property of the pure Einstein gravity, while this can be 
destroyed by introducing complicated couplings that creates Brans-Dicke hairs. 

2. M > 0 limit 

Then unless (3 = 0, are there any hope to recover a Cauchy horizon? By introducing a potential of the 
Brans-Dicke field, we can adjust the Brans-Dicke hair, and this can be helpful to re-create a Cauchy horizon. 

Comparing with FIG. [2} in FIG. [5] we can vary M from 0.1 to 5. This shows that as we increase the effect 
of M, the causal structure is changed and the inner Cauchy horizon appears again as the case of the usual 
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(D = -l, B—-0.5, M—0 



co=-1.4, B=-0.5,M=0 



co—1.4, B=-l, M=0 



v 



FIG. 3: r (left) and $ (right) by varying /? = —0.5 and ft = —1, while we choose M = 0. In addition, for a comparison, 
we also tested not only oj = —1, but also uo = —1.4. Here, the thick black curves are singularities in terms of the Einstein 
frame. 


r vs r along constant V <!V vs U along constant V 




FIG. 4: r (left) and 'Fr 2 (right) for /3 = —0.5 and oj = —1 case. Although r increases to a large value for an in-going 
null observer, the physical areal radius in the Einstein frame goes to zero (hence, towards a singularity). 
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o>=-l,|i=0.5,M=2 cd=- 1, fi=l, M=2 



v V 


(u=-l,P=0.5,M=3 to=-l. 6=1. M=3 



40 ■ ^ V 


-- 
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FIG. 5: r by varying M = 0.1, 1, 2, 3, and 5 (from top to bottom) for /3 = 0.5 (left) and 1 (right). 

charged black holes with Einstein gravity. If we compare with $ (FIG. [6]), then we can see clear dependence. 
Because of the choice of /3, after a formation of a black hole, the Brans-Dicke field is biased toward a weak 
coupling limit, and this is the same as that of the case of M = 0. However, because of the mass parameter, the 
sensitivity of the Brans-Dicke field is limited. Therefore, if the mass parameter is small enough (e.g., M = 0.1) 
so that the Brans-Dicke field is biased to the weak coupling limit even outside the outer apparent horizon, 
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FIG. 6: $ by varying M = 0.1, 1, 2, 3, and 5 (from top to bottom) for j3 = 0.5 (left) and 1 (right). 

then the effects of charge is already screened and hence the internal structure does not have a Cauchy horizon. 
However, as M increases, e.g., for the case of M = 1 or 2, although the Brans-Dicke held inside the horizon is 
biased toward the weak coupling limit, it is no more biased in terms of the outside the outer apparent horizon. 
If we choose much larger value of M (e.g., M = 5), then even the inside the outer apparent horizon, the bias 
of the Brans-Dicke held is limited and decreases. 
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U 


co=-l,|3=1.5,M=3 



v 



u 



FIG. 7: r (left) and $ (right) for /? = 1.5 by varying M = 0.1, 1, 2, and 3 (from top to bottom). 

In addition, as a check, we observe the case when /3 > 1 and varying M in FIG. [7] As we expected, for 
example if (3 = 1.5, during and after the gravitational collapse, the Brans-Dicke field should be biased to weak 
coupling direction, and hence there should be no Cauchy horizon. However, by increasing M , the Brans-Dicke 
hair is controlled and eventually there appears a Cauchy horizon and mass inflation. 


IV. DISCUSSION 

In this paper, we investigated dynamics of string-inspired charged black holes. In our previous paper ;10], 
we already investigated the causal structures and responses of the Brans-Dicke field. Regarding the causal 
structure, the main points were as follows: (1) if /3 > 0, then after the gravitational collapse, there appears a 
Brans-Dicke hair that is biased toward a weak coupling limit and (2) unless /3 = 0, there is no Cauchy horizon 
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and no mass inflation inside a charged black hole. 

In this paper, we investigated the physical origin of the existence or absence of mass inflation and the Cauchy 
horizon. What we have observed are as follows: 

- If there is no potential of the Brans-Dicke field, unless j3 = 0, after the gravitational collapse, the Brans- 
Dicke hair will be formed and biased either toward a weak coupling limit (if /3 > 0) or a strong coupling 
limit (if (3 < 0). For all cases, there is no Cauchy horizon and mass inflation; hence, the existence/absence 
of a Brans-Dicke hair is related to the existence/absence of mass inflation and Cauchy horizon inside a 
charged black hole. 

- Even though /3 ^ 0, by introducing a potential of the Brans-Dicke field, we can reproduce a mass inflation 
and a Cauchy horizon. In other words, if the Brans-Dicke field is well trapped by a potential, then even 
though the Brans-Dicke field is affected by charges of the black hole, we can still see a mass inflation and 
a Cauchy horizon. 

Therefore, we can qualitatively conclude that the existence of mass inflation is related to the existence of the 
Brans-Dicke hair. In other words, we may further say that (within some limited conditions) the existence of 
mass inflation and a Cauchy horizon inside a charged black hole is a consequence of no scalar hair. 

In string theory, the dilaton field should be well trapped by a potential; unless it is unrealistic. Therefore, 
it may be sound to conclude that there may exist mass inflation inside realistic charged black holes in our 
universe. On the other hand, if the trap of the dilaton field is not deep enough and hence a kind of scalar hair 
is detectable by an experiment, then we may conclude that the black hole just have a space-like singularity 
without a Cauchy horizon. 

It may be interesting that we can know the internal information of a black hole by observing the dilaton 
behaviors outside the horizon. However, our descriptions are still qualitative. To say more details, analytic 
modeling and calculations would be required to give a detailed connections between the dilaton hair shape and 
the internal structures. In addition, the dependence or sensitivity by varying another parameters, for examples 
space-dimensions, topology, or the background cosmological constant, can be an interesting future topic. We 
remain them for a future research topic. 
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Appendix A. Implementation to double-null formalism 

We use the double-null coordinates 


ds 2 = —a 2 (u,v)dudv + r 2 (u,v)dfl 2 


(18) 


and present every equations using this metric ansatz. To simplify and present every equations by first order 
differentials, we define the following variables [una: The metric function a, the radial function r, the Brans- 
Dicke field $ and a complex scalar field s = \[4iT<f>, and define 


h = 


a 


7 

d = —, 
a 


f = r, u , g = r }V , W = Z = w = s t . 


z = s, 


(19) 


Using this, the Einstein tensor and the energy-momentum tensor for the Brans-Dicke part and the scalar field 
part are as follows: 


G U u 
G U v 
G vv 

G eS 

pBD _ 


— (/,«-2 fh), 

r 

^2 + « 2 + 4 fa) 

2 

-- ( 9,v ~ 2gd ), 


4 r A 

■~4—^ d u -)- 

cr V r 


1 

87r>I> 
Z 


(W,u - 2 hW) 


p BD __ 


87T$ 2 

gW + fZ a 2 V 


W 2 , 


pBD 


rp BD _ 

1 06 — 


p^ _ 


pKy _ 


pL> _ 


7^c _ 

1 ee — 


87T$ 
1 


47T 


327t$ 


87t4> 


(Z „ - 2 dZ) + 


oo 


87r4> 2 


27ra 2 4> ' u 


47ra 2 ^ 


(gW + fz) 


J^^WZ- — 

47r$ 2 a 2 167t$ ’ 


\ww + iea(ws — ws ) + e 2 a 2 ss] , 


47T 

(a «) 2 

47ra 2 

b**' 


47ra 2 


(wz + zw) + iea(zs — zs) + 


2 (a,, 


ar 


where q(u,v ) = 2 r 2 a jV /a 2 is interpreted as the charge function. 


( 20 ) 

( 21 ) 

( 22 ) 

(23) 

(24) 

(25) 

(26) 

(27) 

(28) 

(29) 

(30) 

(31) 
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After a simple calculation, we can represent equations for a tUV , r tUV , and the field equation for $. We define 
X = <&@X for any quantity X. Then, 


' ( l °S a ),uv 


1 

^,uv 

~ )i2 

y ^ ,uv j 

V 


' ' 2<S> 

0 r 2 
0 0 r 


\ 

' 21 ' 



/ 

U/ 


(32) 


r,w = 2gd - ^(Z, v - 2 dZ) - Z 2 - 


where 


8 $ ’ 


21 - -^oe-l^W + fZ)- ^WZ- a ' 2V 


<£= -fZ-gW- 


2irrcx 2 / ~ c 


3 4- 2w 


- 2 fiC 


'EM 


ra 


4(3 + 2 uj) 


($V" - 2V ). 


In addition, matter Held equations are 

a 2 q 
~TV 2 ' 


a.v = 


ier -x o Z 

Q,v = -~^-(sz - sz) - 0q — , 


&.vv — 


Q,u = 


(“-?) 


ieor cril 

'-zlQ -— ( zs ~ sz ) - Pd 


r J 4 

ier 2 2 2 - o W 

-(sit; — sw) — re ass — pq 

2 J $ 


2 / 

2 r A$’ 


9- 


^. UV 


a 2 A /\ *ea 2 , , a 2 o „ a 2 IT 

r 2 \ r 
/3 gw iearz 
r 


2 , _ a 2 2 _ . 

4 -«*>)- ass-/3 9 - 2r2$ , 

ieags ie 2 /3 . . 

- ; a gs-(1V2 + Zw + zesaZ ). 

r 4r 2 2$ 


r r 

We used numerical integration via a standard 4th order Runge-Kutta method. Consistency 
tests are discussed in Appendix B. 


(33) 

(34) 

(35) 

(36) 

(37) 

(38) 

(39) 

(40) 

(41) 

(42) 

(43) 

and convergence 


Appendix B. Boundary conditions and free parameters 

We specify boundary conditions of all variables (a,r, $,s,a) on the initial u = u\ and v = Vi null surfaces. 
For all simulations in this paper, we choose iq = Vi = 0 and computational domain to have size v = [0; 60] and 
u = [0; 40]. 

By using the gauge freedom to choose r, we choose r(u, 0) iU = r u o < 0 and r(0, v) yV = r v o > 0 such that the 
function r for an in-going observer decreases and that for an out-going observer increases. In addition,we fix 
the effective gravitation constant G = l/4> as unity and we set '!>('«, 0) = 4>(0, v) = 1. 

- In-going null direction : To assign all functions for the in-going null direction, we simply choose s(u, 0) = 0, 
a(u, 0) = 1, q(u,0) = 0, and a(u, 0) = 0, so that the interior is not affected by the gravitational collapse. 
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Since the Misner-Sharpe mass function is 




(44) 

if we give a condition that this vanishes at u = v = 0, it is convenient to choose r iU (w, 0) = —1/2 and 
r iV (0,v) = 1/2. In addition we further choose r(0,0) = ro = 20. In addition, by using the constraint 


equations (Equation (33)), we can complete the assignment of initial conditions. 

- Out-going null direction: We choose a function s(0,u) to induce a collapsing shell. We use 


s(0, v) = A sin 4 ( 7r 


V{ — Vi 


cos 27 t 


Vf - Vi 


V — %)' 

+ i cos ( 2-7T--— ttS 

Vf - Vi 


(45) 


for Vi < v < Vf and otherwise s(0,u) = 0. By using this, we integrate Equation (34) to determine a(0,i>) 


on the u = 0 surface. Also, at the same time, we integrat Equations (38) and (39) to determine g(0,u) 
and a(0, v). 


Now for a pulse, we choose Vf = 20, A = 0.15, and 6 = 0.5, leaving four free parameters (uj, /3, e, M), where uj 
is the Brans-Dicke coupling parameter, /? is the coupling between the matter sector and the Brans-Dicke field, 
e is the gauge coupling, and M is the mass parameter of the potential of the Brans-Dicke field. Note that for 
each choice of model parameters, we assign meanings as follows m- 


— uo : 0 ( f(R ) limit), —1 (dilaton limit), —1.4 (braneworld limit), —1.6 (ghost limit), 


— /3 : 0 (Type HA), 0.5 (Type I), 1 (Heterotic), 

— e : 0 (neutral black hole), 0.3 (charged black hole). 


Overall this paper, we fix e = 0.3 (only consider charged cases), and vary uj, /3, and M. 

We solved this system by using the same code with the authors’ previous paper m- We repeat a demon¬ 
stration of convergence and consistency checks: FIG. [8| Here, we define 


£ ( x n) = 


U N 


b 2N I 


X 


HighRes I 


(46) 


where x l N is the i-th grid of a dynamic variable x with resolution N; xJjjgiiRes numerical 


resolution what we used. In addition, to check constraint equations (Eqs. (33) and (34)), we define 

\Ck\ 


X (Cat) = 


\Gh 


(47) 


HighRes I 


where C' l N is the residual of constraint equations (let us say that C uu or C vv ) at the *-th point for simulation 
with resolution N and GjjighRes a corresponding Einstein-tensor component ( G uu or G vv ). We can see good 
convergence and consistency from FIG. |8} as resolution increases, errors are well controlled and converges less 
than 0.1 %. 
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FIG. 8: Convergence and consistency checks for along a line u = 20 for a simulation with parameters: to = —1.4, e = 0.3, 
p = —1.0, and M = 0. As the resolution increases, errors are well controlled less than 0.1 %. 


Appendix C. Catalog of metric function and energy-momentum tensors 

In this appendix, we summarize numerical data of a 2 , T uu , and T vv . The conditions that we have used are 
as follows: 


- String limit, where we use /? = 0.5 and /? = 1, u> = —1 (dilaton limit), and varying the mass parameter 
M = 0.1, 1, 2, 3, and 5: FIGs.[9j[T0j and [IT 
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- f) < 0, where we use /3 = —0.5 and —1, by fixing M = 0 and varying w = — 1 and —1.4: FIGs. 12 13 and 


El 


- (3 > 1, where we use (3 = 1.5, by fixing ui 


— 1 and varying M = 0.1, 1, 2, and 3: FIGs. 15 and 16 


From these data, when there is a Cauchy horizon, we can see that a approaches to zero, where this is an 
evidence of the existence of mass inflation, as was observed by |15| (for more justification, see also Appendix A 

of HD]). 
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FIG. 13: T uu by varying /? = —0.5 (upper), —1 (lower) and varying u> = —1 (left), —1.4 (right). 
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